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Abstract 

Eulcr's celebrated buckling formula gives the critical load N for 
the buckling of a slender cylindrical column with radius B and length 
L as 

iV/(7r3B2) = (£/4)(i?/L)2, 

where E is Young's modulus. Its derivation relies on the assumptions 
that linear elasticity applies to this problem, and that the slenderness 
(B/L) is an infinitesimal quantity. Here wc ask the following question: 
What is the first nonlinear correction in the right hand-side of this 
equation when terms up to (B/L)^ are kept? To answer this question, 
we specialize the exact solution of incremental non-linear elasticity for 
the homogeneous compression of a thick compressible cylinder with 
lubricated ends to the theory of third-order elasticity. In particular, 
we highlight the way second- and third-order constants — including 
Poisson's ratio — all appear in the coefficient of (B/L)^. 
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1 Introduction 

One of the first, and most important, problems to be tackled by the theory of 
linear elasticity is that of the buckling of a column under an axial load. Using 
Bernoulli's beam equations, Euler found the critical load of compression A^^cr 
leading to the buckling of a slender cylindrical column of radius B and length 
L. As recalled by Timoshenko and Gere [T], Euler looked at the case of 
an ideal column with (built in)-(free) end conditions. What is now called 
"Euler buckling", or the ^^fundamental case of buckling of a prismatic bar" 
[1] corresponds to the case of a bar with hinged-hinged end conditions. The 
corresponding critical load is given by 



where E is the Young modulus. The extension of this formula to the case of a 
thick column is a non-trivial, even sophisticated, problem of non-linear three- 
dimensional elasticity. In general, progress can be only made by using reduc- 
tive (rod, shells, etc.) theories. However, there is another choice of boundary 
conditions for which the criterion ([T]) is valid: namely, the case where both 
ends are "guided" or "sliding" (the difference between the two cases lies in 
the shape of the deflected bar, which is according to the half-period of a sine 
in one case and of a cosine in the other case.) In exact non-linear elastic- 
ity, there exists a remarkable three-dimensional analytical solution to this 
problem (due to Wilkes [3J) which describes a small-amplitude (incremen- 
tal) deflection superimposed upon the large homogeneous deformation of a 
cylinder compressed between two lubricated platens. In this case, the Euler 
formula can be extended to the case of a column with finite dimensions, for 
arbitrary constitutive law. 

The exact incremental solution allows for an explicit derivation of Euler's 
formula at the onset of non-linearity, which combines third-order elastic con- 
stants with a term in (B/L)^. In Goriely et al. [2], we showed that for an 
incompressible cylinder. 



where A is Landau's third-order elasticity constant. This formula clearly 
shows that geometrical non-linearities (term in (B/L)'^) are intrinsically cou- 
pled to physical non-linearities (term in A) for this problem. (For connections 
between Euler's theory of buckling and incremental incompressible nonlinear 
elasticity, see the early works of Wilkes [3j , Biot \^ , Fosdick and Shield [S] , 
and the references collected in [2].) 
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Now, in third-order m compressible elasticity, there are two elastic con- 
stants [6]: the shear modulus fi {= E/3) and A. In third-order compressible 
elasticity, there are five elastic constants: A and fi, the (second-order) Lame 
constants (or equivalently, E and u, Young's modulus and Poisson's ratio), 
and A, B, and C, the (third-order) Landau constants. Euler's formula at 
order {B/Ly, equation ([T]), involves only one elastic constant, E. It is thus 
natural to ask whether Poisson's ratio, z/, plays a role in the non-linear cor- 
rection to Euler formula of (B/L)^, the next-order term. The final answer 
is found as formula (32) below, which shows that the non-linear correction 
involves all five elastic constants. 



2 Finite compression and buckling 

In this section, we recall the equations governing the homogeneous compres- 
sion of a cylinder in the theory of exact (finite) elasticity. We also recall the 
form of some incremental solutions that is, of some small-amplitude static 
deformations which may be superimposed upon the large compression, and 
which indicate the onset of instability for the cylinder. 

2.1 Large deformation 

We take a cylinder made of a hyperelastic, compressible, isotropic solid with 
strain-energy density W say, with radius B and length L in its undeformed 
configuration. We denote by [R, 0, Z) the coordinates of a particle in the 
cylinder at rest, in a cylindrical system. 

Then we consider that the cylinder is subject to the following deformation, 

r = Aii?, ^ = e, z = hZ, (3) 

where (r, 6, z) are the current coordinates of the particle initially at (i?, G, 
Z), Ai is the radial stretch ratio and A3 is the axial stretch ratio. Explicitly, 
Ai = h/B and A3 = Z/L, where h and / are the radius and length of the 
deformed cylinder, respectively. The physical components of F, the corre- 
sponding deformation gradient, are: F = Diag (Ai, Ai, A3), showing that the 
deformation is equi-biaxial and homogeneous (and thus, universal). 

The (constant) Cauchy stresses required to maintain the large homoge- 
neous compression are (see, for instance, Ogden [7]), 

or, = (det F)"^ XiWi, i = 1, 2, 3 (no sum), (4) 

where Wi = dW/dXi. In our case, cxi = (T2 because the deformation is equi- 
biaxial, and 0"! = a2 = because the outer face of the cylinder is free of 
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traction. Hence 



ai=a2 = Ar'Ag Wi = 0, = X['W^. (5) 

Note that we may use the first equahty to express one principal stretch in 
terms of the other (provided, of course, that inverses can be performed). 



2.2 Incremental equations 

Now we recall the equations governing the equilibrium of incremental solu- 
tions, in the neighborhood of the finite compression. They read in general as 

m 

div s = 0, (6) 

where s is the incremental nominal stress tensor. It is related to u, the 
incremental mechanical displacement, through the incremental constitutive 
law, 

s = B (grad u)^ , (7) 

where B is the fourth-order tensor of incremental elastic moduli and the 
gradient is computed in the current cylindrical coordinates [8] . The non-zero 
components of B, in a coordinate system aligned with the principal axes 
associated with the deformation ([s]), are given in general by [7] 

i 7^ J, ^ Xj, 

i 7^ J, \ 7^ 

i ^ j, \i = Xj, 
i ^ j, Xi = Xj, (8) 

(no sums), where Wij = d'^W/{dXidXj). Note that here, some of these 
components are not independent one from another because Ai = A2 and 
o"i = o"2 = 0. In particular, we find that 

-B1212 = -B2121 = -B12215 -S2323 = -S1313 = -B1331 = -S2332, 

-B2222 = -Bull, -B2233 = -Bll33) -B3232 = -B3131, -B1122 = -Biiii — 2i?i2l2- 

(9) 



JBiijj — XiXjWij, 



JB,,,, = {Xm - X,W,)X^/iX^ - X]), 
JBijj, = {X,W, - X,W,)X,X,/{Xl - A|), 

J -Bijij i^Biiii -Biijj ~\~ ^i^^'i) /'^•) 



2.3 Incremental solutions 

We look for incremental static solutions that are periodic along the circum- 
ferential and axial directions, and have yet unknown radial variations. Thus 
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our ansatz for the components of the mechanical displacement is the same 
as Wilkes's [3]: 

Ur = Ur{r) cos n9 cos kz, uq = Uelr) sin n6 cos kz, = Uz{r) cos nO sin kz, 

(10) 

where n = 0, 1, 2, . . . is the circumferential mode number; k is the axial 
wavenumber; the subscripts {r,6,z) refer to components within the cylin- 
drical coordinates {r,6,z); and all upper-case functions are functions of r 
alone. 

Dorfmann and Haughton [S] show that the following displacements U^^\ 
U^^\ and U^^^ are solutions to the incremental equations ([6]), 



and 



Il,{qkr) 



■n J , , . (filing^ - B: 
In{qkr), 



3131; 



qkr 



q{Bisi3 + -B1133) 



In{qkr) 



-Iniqskr), -—I'^{qskr), 
r n 



(11) 
(12) 



where q = qi,q2 and In is the modified Bessel function of order n. Here qi, 
q2, and q^ are the square roots of the roots qf, g| of the following quadratic 
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m q 

-Bi3i3i?iiiig^+ [(-B1133 + -81313)^ — -Bi3i3i?3i3i — -B3333-Bllll]q'^ + -B3333-B3131 = 0, 



and of the root of the following linear equation in 

Bi2i2q^ — -B3131 = 0, 



(13) 
(14) 



respectively. 

From ([7]) we find that the incremental traction on planes normal to the 
axial direction has components of the same form as that of the displacements, 
namely 



Srr = Srr{r) COSTlO COskz^ 

Sr-e = Sr-eir) sinnO cos kz, 
Sry = SrAr) cosuOsinkz, 



(15) 



say, where again all upper-case functions are functions of r alone. Then we 
find that the traction solutions corresponding to the solutions (11)-(12) are 



5 



given by 

rS'"(r), rS'^>(r) = {2Bui2lL{<lkr) 



V Qkr g (5i3i3 + 5ii33) y 

2nSi2i2 ^ LiQkr) , -M1313 — — + 1 4 g/cr 

\ qkr J \ Bi3i3 + ^1133 / 



T 



(16) 



and 



2i?i2i2 ( ^^^^ - qskl'^iq^kr) 



\ n \ r n J J 

(17) 

The general solution to the incremental equations of equilibrium is thus 
of the form 

r5(r) = [r5W(r) | rS^^\r) \ rS^'^\r)\ c, (18) 

where S = [Srr, Sre, SrzY' c is a constant three-component vector. Note 
that we use the quantity rS for the traction (instead of S), because it is the 
Hamiltonian conjugate to the displacement in cylindrical coordinates [9]. 

Now when the cylinder is compressed (by platens say), its end faces should 
stay in full contact with the platens so that the first incremental boundary 
condition is 

u, = 0, on z = 0,l, (19) 

which leads to [HI, |2] 

k = mvr//, (20) 

for some integer m, the axial mode number. From (15), we now see that on 
the thrust faces, we have 

Srz = 0, on z = 0,l, (21) 

which means that the end faces of the column are in sliding contact with 
the thrusting platens. In other words, in the limit of a slender column, we 
recover the Euler strut with sliding-sliding, or guided-guided end conditions. 
In Figure [T| we show the first two axi-symmetric and two asymmetric modes 
of incremental buckling. 

The other boundary condition is that the cylindrical face is free of incre- 
mental traction: S{h) = 0. This gives 

A = det [bS^^\b) I 65^') (6) | bS^^\b)] = 0. (22) 
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Figure 1: First two axi-symmetric and two asymmetric modes of buckling 
for a compressed strut with guided-guided end conditions: n is the circum- 
ferential mode number and m the axial mode number. For slender enough 
cylinders, the n = 1, m = 1 mode is the first mode of buckling. 

3 Euler buckling 

3.1 Asymptotic expansions 

We now specialize the analysis to the asymmetric buckling mode n = 1, m = 
1, corresponding to the Euler buckling with guided-guided end conditions, 
in the limit where the axial compressive stretch A3 is close to 1 (the other 
modes are not reached for slender enough cylinders). To this end, we only 
need to consider the so-called third-order elasticity expansion of the strain 
energy density, for example that of Landau and Lifshitz [10] , 

W = ^ (tiEf + ft tr(E^) + ^tr(E3) + B (trE) ti{E^) + ^ (trE)=^ , (23) 
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where E = E"^ is the Lagrange, or Green, strain tensor defined as E = 
(F"^F — I) /2, A and yU are the Lame moduh, and A, B, C are the Landau 
third-order elastic constants (Note that there are other, equivalent, expan- 
sions based on other invariants, such as the ones proposed by Murnaghan 
Toupin and Bernstein |12j. Bland [13j, or Eringen and Suhubi [lAl, see 
Norris [15] for the connections.) 

To measure how close A3 is to 1, we introduce e, a small parameter pro- 
portional to the slenderness of the deformed cylinder, 

e = kb = nb/l. (24) 

Then we expand the radial stretch Ai and the critical buckling stretch A3 in 
terms of e up to order M, 

M M 

Ai = Ai(e) = l+^ape^+(9(e^^+i), A3 = A3(e) = 1 + ^ /3pe^+0(e^^+i), 
p=i p=i 

(25) 

say, where the a's and /3's are to be determined shortly. Similarly, we expand 
A in powers of e. 

Mi 

p=i 

say, and solve each order c/p = for the coefficients Up and /3p, making use of 
the condition ai = 0. We find that Up and (5p vanish identically for all odd 
values of p, and that Ai and A3, up to the fourth-order in e, are given by 

Ai = l + «2e' + «4e^ + 0(e^), \^ = I + ^e^ + (3^6^ + 0{e^), (26) 

with a2 and 0:4 given by 

V 

v{\ + v) 

= 

32 

_ (1 + //)(! -2;/) 
16E 

wherein 

^' = -1 

_29 + 39z/ + 8z/2 

P4 - 



[u'^A + (1 - 2z/ + 6z/2)i3 + (1 - 2ufC] - (27) 



96(1 + 1/) 

+ [(1 - 2iy^)A + 3(1 - 2z/)(l + 2z/2)i3 + (1 - 2iyfC] . {2t 
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Note that we switched from Lame constants to Poisson's ratio and Young's 
modulus for these expressions, using the connections v = A/(2A + 2fi) and 
E = /i(3A + 2/i)/(A + /i). 



3.2 Onset of nonlinear Euler buckling 

The analytical results presented above are formulated in terms of the current 
geometrical parameter e, defined in (24). In order to relate these results 
to the classical form of Euler buckling, we introduce the initial geometric 
slenderness B/L. Recalling that e = nb/l, A3 = l/L, and b = XiB, we find 
that 

eA3 = vrAi(5/L). (29) 



We expand e in powers of B/L, and solve (29) to obtain 

e = n{B/L) + {a2- P2)TT\B/L)^ + 0{iB/L)^) 
= n{B/L) + (1 + u){ny4){B/Lr + O {{B/Lf) . 



(30) 



Second, we wish to relate the axial compression to the current axial load 
A^. To do so, we integrate the axial stress over the faces of the cylinder. 



= -27r 



ra^dr 



-nXlB^as, 



(31) 



because is constant, given by ([5])2. 

Finally, in order to write the nonlinear buckling formula, we expand Ai 
and A3 in (31), using (26), and then expand e in powers of the slenderness 
{B/L), using (30). It gives the desired expression for the first non-linear 
correction to the Euler formula. 



TT^B^ 



E 



96~ 
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NL 



(32) 



where 



NL 



, 13 + 12i/ - 2z/2 



E 



+ 12 [(1 - 2u^)A + 3(1 - 2z/)(l + 2u^)B + (1 - 2u)^C] . (33) 

We now check this equation against its incompressible counterpart 
Theoretical considerations and experimental measurements [6], [T7], [18], 
|19j . show that in the incompressible limit, E and A remain finite, u — )• 1/2, 
(1 — 2i>)B — )• —E/3, and (1 — 2z/)^C — )■ 0. It is then a simple exercise to verify 
that (32) is indeed consistent with ^ in those limits. 
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3.3 Examples 



Table 1: Lame constants and Landau third-order elastic moduli for five solids 



(10^ N- m-^) 



material 


A 




A 


B 


C 


Polystyrene 


1.71 


0.95 


-10 


-8.3 


-10.6 


Steel Hecla 37 


111 


82.1 


-358 


-282 


-177 


Aluminium 2S 


57 


27.6 


-228 


-197 


-102 


Pyrex glass 


13.5 


27.5 


420 


-118 


132 


Si02 melted 


15.9 


31.3 


-44 


93 


36 



To evaluate the importance of the non-linear correction, we computed the 
critical axial stretch ratio of column buckling for two solids. In Table [l| we 
list the second- and third-order elastic constants of five compressible solids, 
as collected by Porubov [16J (in the Table we converted the "Murnaghan 
constants" given by Porubov to the Landau constants A,B,C). Figure [2] 
shows the variations of A3 with the squared slenderness (B/L)"^, for pyrex 
and silica (two last lines of Table 1). 

4 Conclusions 

The present analysis provides an asymptotic formula for the critical value of 
the load for the Euler buckling problem, with guided-guided (sliding-sliding) 
end conditions. This formula was checked both in the incompressible limit 
and on particular cases against the exact value of the buckling obtained from 
the exact solutions. Not surprisingly it reinforces the universal and generic 
nature of the Euler buckling formula as the correction is small for most 
systems even when nonlinear elastic effects and nonlinear geometric effects 
are taken into account. It would be of great interest to see if these effects 
could be observed experimentally. 
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